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Abstract: We study a number of U(l)x models featuring a Dirac fermion dark matter 
particle. We perform a comprehensive analysis which includes the study of corrections to 
the muon magnetic moment, dilepton searches with LHC data, as well as direct and indirect 
dark matter detection constraints. We consider four different coupling structures, namely 
U(l) B -L,U(l) d - u ,U(l) universa i, and t/(l) 10+ g, all motivated by compelling extensions to 
the standard model. We outline the viable and excluded regions of parameter space using 
a large set of probes. Our key findings are that (i) the combination of direct detection 
and collider constraints rule out dark matter particle masses lighter than ~ 1 TeV, unless 
rather suppressed Z'-fermion couplings exist, and that (ii) for several of the models under 
consideration, collider constraints rule out Z' masses up to ~ 3 TeV. Lastly, we show that 
we can accommodate the recent Diboson excess reported by ATLAS collaboration within 
the U(l) d -u model. 
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1 Introduction 

The existence of dark matter (DM) has been ascertained by numerous lines of evidence, 
originating from different scales and epochs in the history of the universe. However, the 
fundamental particle nature of DM is yet to be unveiled. From a particle physics standpoint, 
thermally produced DM can have have mass in the keV-TeV range. Within this mass 
range, WIMPs (Weakly Interacting Massive Particles) are some of the most compelling 
dark matter candidates (see for example the reviews in Ref. [1-7]). 

WIMPs have been intensively searched for in experiments that use a broad variety of 
techniques, customarily grouped in the three categories of (i) direct detection, (ii) indirect 
detection, and (iii) collider searches [8]. 

(i) Direct detection relies on the measurement of the energy recoil deposited by WIMPs 
on nuclei in underground experiments typically operating at very cold temperatures. With 
a good control over the background, current experiments seek to observe excesses of scat¬ 
tering events, which are potentially translated to measurements of the WIMP-nucleon scat¬ 
tering cross section and mass, for a given velocity distribution. The difficulty in precisely 
assessing the number of background events has in the past misled several direct detec¬ 
tion experiments to claim the observation of excesses which could be plausibly explained 
by 7 — 30 GeV WIMPs [9, 10]. However, null results from other experiments, including 
CDMSlite, XENON, and LUX [11, 12] are in blatant tension with this possibility, despite 
theoretical efforts to accommodate a WIMP interpretation by means of isospin violation 
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[13-15] or scattering with impurities in the detectors [16]. At present, no conclusive signal 
from WIMP dark matter has been reported in direct detection experiments, and increas¬ 
ingly strong limits on the nucleon-WIMP cross section probe theories beyond the Standard 
Model with unprecedented accuracy. 

(ii) Indirect detection refers to the observation of DM annihilation products, such as 
cosmic- and gamma-rays (see Ref. [17] for a pedagogical introduction). This is an exciting 
and promising avenue, which makes use of a variety of different messenger particles and 
of observational devices, and which relies on direct or indirect information about the DM 
density distribution in the Galaxy and beyond. Ref. [18] provides an overview of recent 
results. Since potential discoveries in indirect detection generically suffer both from large 
and vastly unknown background and systematic uncertainties, we focus here on bounds 
derived using the relatively robust limits currently set by Fermi-LAT observations of dwarfs 
spheroidal galaxies [19]. 

(iii) Collider searches for dark matter are often based on the search for events with 
large missing energy, associated with pair-produced WIMPs escaping the detector [20-22], 
In models where the mediator between the visible and WIMP sector has sizable couplings 
to SM fermions, either dijet or dilepton limits are typically most stringent [23-52], whereas 
mono-jet ones are complementary, and often competitive in the low WIMP mass regime 
[53-58]. As explained below in Sec. 5, the dilepton searches provide the most stringent 
collider bounds for the models that we consider in this study. 

Of additional relevance for the class of models under discussion here are limits from 
contributions to the muon magnetic moment, especially in the regime where the DM cou¬ 
plings to the mediator are suppressed, i.e. in a regime where direct detection and indirect 
detection limits weaken. Here, we use an adapted version of the public code [59] to assess 
such limits. 

In view of the current bounds, several DM portals have been studied with the purpose 
of outlining the remaining viable parameter space after combining existing limits. Some of 
the most interesting and studied candidates are Dirac fermions. In this context, Z' gauge 
bosons are natural mediators since they appear in a large multitude of gauge extensions 
of the SM. Most often, this Z' portal is studied in simplified models such as Ref. [24— 
26]. The simplified model approach is valid and interesting. However, the conclusions 
drawn often cannot be directly applied to particular models, due to the large variations 
in the structure and magnitude of the vector and axial-vector couplings. Therefore, it is 
worthwhile to study particular gauge structures, which appear in several extensions of the 
SM. These types of extensions often address intriguing open questions such as neutrino 
masses, number of family generations, and matter and anti-matter asymmetry. 

An additional motivation to focus on Z' models originates from the recent excess 
reported by the ATLAS collaboration in the resonant diboson to hadronic final states 
channel at an energy of around 2 TeV [60]; such excess, currently featuring a local statistical 
significance of 3.4er, 2.6cr and 2.9er respectively in the WZ, WW and ZZ channels [60], has 
been ascribed to an extra 2 TeV bosonic particle [61-63]. In particular, Ref. [62] constructs 
an explicit leptophobic 17(1)' extension inspired by an Eq supersymmetric GUT. While we 
do not consider this specific model here, our results provide a direct connection between the 
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ATLAS diboson excess and models with a dark matter particle candidate, and a roadmap 
for further studies in this direction. 

In summary, the current experimental landscape allows us to probe models of new 
physics with a WIMP dark matter candidate from several independent directions. Here, 
for the first time, we perform a comprehensive analysis taking advantage of collider, (g— 2)^, 
direct, and indirect WIMP dark matter detection limits for four theoretically simple and 
well-motivated U (1) gauge extensions. In particular, we outline the remaining viable region 
of parameter space and that which is now experimentally excluded, for a broad range of 
WIMP masses ranging from 8 GeV up to 5 TeV. 

The remainder of this study is as follows: in the following section we describe and 
define the four 17(1) gauge extensions that we consider throughout our study; in Sec. 3 we 
outline the dark matter phenomenology (relic density, direct, and indirect searches) in these 
scenarios; Sec. 4 explores the impact of contributions to the muon anomalous magnetic 
moment, while Sec. 5 estimates the impact of collider searches. Section 6 summarizes and 
compiles all of our results, while Sec. 8 concludes. 

2 17(1) Extensions 

One of the most minimal ways of extending the SM gauge structure is to include an 
additional spontaneously broken U(l)x gauge group. As a result, a new neutral, massive 
spin-1 gauge boson ( Z') arises. In the scenarios of interest here, the DM is a new Dirac 
fermion charged under the U(l)x group. As far as the setup’s phenomenology is concerned, 
the relevant parameters of the theory are the Z' mass (M^/), its coupling to fermions 
(including the DM), its width (Tz')-, as well as the DM mass (M x ). 1 In this work we focus 
on extensions of the SM which (i) avoid flavor changing neutral currents at tree level, and 

(ii) allow for cancellation of triangle anomalies by the introduction of vector-like fermions. 
It has been shown that these two requirements greatly restrict the list of possible coupling 
structures. We will focus on four interesting cases [67, 68], namely: 

(i) 1/(1 )b~l'- baryon and lepton numbers are accidental anomalous global symmetries 
of the SM, whereas B-L is not. Thus an interesting extension of the SM consists of gauging 
the B-L symmetry [67]; 

(ii) U(l)d- U - this structure arises after spontaneous symmetry breaking in some E$ 
grand unified theories [69]; 

(iii) L/(l) 1 o-t- 5 : this structure arises after the spontaneous symmetry breaking of flipped 
SU(5) grand unified theories [68]. 

(iv) U(l) un i versa i: an extensively explored extension which inherits the universal cou¬ 
pling structure of the SM Z [67]. 

As mentioned above, we seek suppression of flavor-changing neutral current processes 
at tree level, thus we focus on models where the Z' couplings to the fermions are generation- 
independent [68]. The relevant charges of the left-handed quark doublet (Ql), right-handed 
quarks (), lepton doublet II, and right-handed leptons Ir are shown in Table 1. 

1 Some recent model-independent attempts to study the DM phenomenology in the context of U(l)x 
theories include Ref. [24-26, 36, 39, 64-66]. 
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Table 1 . The generation-independent exotic Z' 17(1) charge assignments explored in this paper. In 
all cases there is a free continuous parameter x which, for simplicity, we take equal to 1 throughout 
our analysis. 


In what follows, we do not assume any specific scalar structure inducing spontaneous 
symmetry breaking of the U{l)x gauge group and that the SM Higgs is neutral under U (l)x 
so that there is no Z — Z' mixing at tree level, which relaxes otherwise strong electroweak 
precision constraints. Finally, since the DM particle is postulated to be charged under 
U(l)x, it will affect the anomaly cancellation requirement; however, in principle additional 
charged fermions might be postulated, with no impact on the DM phenomenology (see 
Ref. [45] for an explicit example). Therefore, with the exception of the Z' and of the DM 
particle, we assume any additional fields beyond the SM is effectively decoupled, and will 
not discuss them further in association with the models’ phenomenology. 

The four U(l)x charge assignments shown in Table 1 depend on a continuous free 
parameter denoted by “x”. Although the successful cancellation of triangle anomalies is 
guaranteed for any value of x, for simplicity, we fix x = 1. We will also restrict the new 
17(1) Z l gauge coupling (gz 1 ) to two numerical values: gz> = 1 or 0.5. With x and gz 1 
specified, we can then determine the numerical values of the Z '-SM fermion couplings as 
discussed in Sec. 3 below. 

3 Dark Matter Phenomenology 

A new massive and neutral Z' is an interesting feature of the low-energy limit of many 
extensions to the SM possessing a viable DM candidate [35, 37, 65, 70-76]. As discussed 
above, the number of possible simple 17(1) charge extensions is greatly reduced after con¬ 
sidering anomaly cancellations and electroweak precision measurements. In particular, we 
focus on the the 17(1) charge assignments in Table 1. We will parametrize the simplified 
Lagrangian responsible for a Z' interacting with a Dirac fermion DM (y) and SM fermion 
(/) as 


^2' [X7 M (ff*v + SW 5 ) X+ Y! ( 3/v + 9fal b ) /] , (3-1) 

/e SM 

where the sum is over all quarks and leptons of the SM (including neutrinos when relevant). 
The Z'-SM couplings above can be specified in terms of the 17(1 )x charges and gauge 
coupling (gz 1 ) since 
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Figure 1. Relevant tree-level processes for: (a,b) pair-annihilation; (c) elastic scattering off of SM 
fermions; (d) muon anomalous magnetic moment. 
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where the z' s are the U(l)x charges of the the gauge eigenstates of Table 1 and gz> is the 
new U( 1) gauge coupling constant. 

We note that our restriction to Dirac DM allows for more general interactions with 
the Z' than Majorana DM, since Dirac DM allows for vector interactions with the Z'. 
Therefore, restricting to Majorana DM instead would affect mostly the scattering behavior 
of DM with quarks inside nuclei since vector interactions are coherent over the entire nucleus 
as will be discussed in Sec. 3.1. We also assume that the y’s vector and axial interactions 
are of the same magnitude and sign, i.e. g xv = g xa . Hence, we define the coupling g x 
as g x = g xv = g xa . The case where the vector and axial couplings have opposite signs, 
g xv = —g X a, would not change any qualitative aspect of our results (see Ref. [26] for studies 
which assume g xv ^ g xa ). 


3.1 Direct Detection 

The general parametrization of the Z’ interactions in Eq. (3.1) induces both spin-independent 
(SI) and spin-dependent (SD) scattering with nuclei. In particular, elastic scattering occurs 
through the t-channel exchange of a Z’ (see Fig.l, c). As mentioned previously, in the case 
where both y and the valence quarks of nucleons possess vector interactions with Z', large 
coherent spin-independent scattering may occur, and this process is severely constrained 
by current bounds from direct detection experiments [11, 77]. An approximate form for 
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the SI cross section at low-momentum transfer is 
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(3.3) 


where /r X n is the WIMP-nucleon reduced mass, and where Z and A are the atomic number 
and atomic mass of the target nucleus, respectively. In scanning through the parameter 
space of different models, we demand that the scattering cross section of Eq. (3.3) be below 
the most stringent current upper bound, given at present by the LUX experiment [11]. 

Similarly, if the nucleon valence quarks and % both possess axial-vector interactions 
with the Z ', one additionally has spin-dependent scattering, which, again in the low mo¬ 
mentum transfer limit, has a cross section of the form 


£j_SD (per neutron) 


3M xneu t 


ala 


7T 


M% V 


g U aK eUt +g d a (Ar + Ar‘) 


1 2 


(3.4) 


where A“ eut are the quark spin fractions of the neutron. We will take these to be A“ eut = 
—0.42, A[) eut = 0.85, A" eut = —0.08 [78]. We then require that the spin-dependent rate to 
be below current published limits from XENONIOO [77]. For spin-dependent bounds, we 
focus on scattering with neutrons since these limits are at present the most stringent. 


3.2 Indirect Detection 

From the Lagrangian of Eq. (3.1), if x possesses non-negligible interactions with Z', then, 
during the early universe, x was i R thermal equilibrium with SM particles. Residual pair 
annihilation can still occur in the late universe, producing SM quark and lepton pairs. In 
turn, this leads to high energy gamma-rays via neutral pion production and final state 
radiation, detectable with gamma-ray telescopes such as the Fermi Large Area Telescope 
[79]. In particular, constraints can be derived from the non-observation of signals in DM- 
dominated targets such as dwarf spheroidal galaxies [19, 80, 81]. 

From the Lagrangian interaction terms indicated in Eq. (3.1), x pairs can annihilate 
through an s-channel Z' exchange process to a pair of SM fermions, as long as M x > mf 
(see (a) of Fig. 1). The non-relativistic form for this annihilation cross section is 


n, 
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(3.5) 


where v is the relative velocity of the annihilating DM pair and n c is the number of colors 
of the final state SM fermion. Sufficiently near resonance, the width of Z', T Z ', should be 
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included in Eq. (3.5). Using again the interaction terms of Eq. (3.1), the Z' width takes 
the simple form 
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where 0 is the unit step function. 

Furthermore, if M x > M^/, then x may also annihilate directly into pairs of on-shell 
Z' bosons, which subsequently decay to SM fermions. The non-relativistic form for this 
annihilation channel is 
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To calculate the gamma-ray yield for this annihilation mode, one needs to account for the 
Z' decay channels. 

Whenever appropriate, we utilize limits on the gamma-ray flux from Fermi observa¬ 
tions. For a given annihilation final state, we thus need to calculate the differential or 
integrated gamma-ray yield. To do this, we utilize the PPPC4DMID code, which calculates 
the spectrum of gamma-rays for direct annihilations to SM fermions [83]. Compared to the 
direct annihilation case, annihilations to pairs of Z' bosons lead to a spread in the gamma- 
ray energy due to the Lorentz boost between the x and Z l rest frames. In this case, the 
spectra from PPPC4DMID, which corresponds to direct annihilations to SM fermions, are 
convolved over a finite energy range, corresponding to the kinematics of the boosted final 
state fermions [84], 

Astrophysical observations place upper bounds on the quantity 
where (av )/ and N~ f are the annihilation rate and number of photons, within a given 
energy-bin, produced in a single annihilation to some final state //, and the factor of 
| takes into account that Dirac DM is not self-conjugate. There are several possible 
annihilation channels whose relative importance depends on the dark matter mass and 
the U (1).y model. In order to derive the indirect detection limits in a consistent way, 
we implement the recent Dwarf PASS8 results from FERMI-LAT [19]. In particular, we 
utilize the upper limits that FERMI Dwarf observations place on direct annihilations to 
bb and rescale accordingly to account for the total photon yield in the Z l model under 
consideration. We point out that competitive but less robust indirect detection limits from 
gamma-ray observations of the center of the Milky Way could be used [82], For other 
complementary limits applicable to this model see Refs. [85-90] 

It is important to note that when M z / <C M x non-perturbative effects such as Som- 
merfeld enhancement might become relevant, especially after freeze-out, and would thus 
significantly strengthen the bounds from indirect detection [91, 92], Since most of the vi¬ 
able parameter space in our models corresponds to larger Z' masses, we do not incorporate 
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this effect into our analysis. We also point out that if had included couplings between the 
Z l and gauge bosons, the shape of the abundance curve would change and mild 7 -ray line 
bounds would be applicable [93]. 


4 The Muon Anomalous Magnetic Moment 

Significant contributions to the muon anomalous magnetic moment a^ are expected in 
U(l)x models where the neutral Z' gauge boson has sizable couplings to leptons. The 
E821 experiment at Brookhaven National Laboratory, which measured the precession of 
muons and anti-muons in a constant external magnetic field as they circulated in a confining 
storage ring, reported the value o ^ 821 = (116592080 ± 63) x 10 -11 [94, 95]. Thus, 


Aa A1 (E821 - SM) = (295 ± 81) x 10“ n , (4.1) 

i.e. an excess of about 3.6 o compared to the SM expectation. This excess might or not 
be associated with new physics, but the result stands in any case as a general constraint 
to models producing significant contributions to a^. 

It is important to notice that due to large uncertainties in the SM hadronic contribu¬ 
tions to a^, caution should be used in taking this excess at face value. The contribution to 
a^ stemming from Z' gauge bosons is of the form [59], 


A a^(Z') 


f 1 , 9%Py(x) + ftgPqjx) 
87T 2 M|, 7 0 T (1 — x)(l — \ 2 x) + A 2 x’ 


(4.2) 


where g and g llA are the vector and axial couplings to muons, respectively, A = m^/Mz' 
and, 


P v (x) = 2x 2 (1 — x) 

P a {x ) = 2x(l — x) ■ {x — 4) — 4A 2 • x 3 . (4.3) 

In the limit Mz> 3> we find [96, 97], 

A aJZ') = ( 1 g 2 uv - 5 g 2 ua ) . (4.4) 

’ 47r 2 M 2 , V3 ^ 3 y/i0 7 v ; 

It is worth pointing out that the overall Z l correction to the muon magnetic moment can 
take positive or negative values depending on the relative magnitude of the vector and 
axial couplings. Rather light Z' gauge bosons are required to explain the excess, and as 
we shall see further those masses are already excluded by collider searches. Nevertheless, 
one can still use the measurement to place lcr limits by demanding the contributions to 
lie within the error bar. This is precisely what we do here. In Table 3, we summarize our 
limits for the four models discussed in our analysis. Those exclusion bounds are also shown 
as horizontal gray lines in Figs. 4-27. We also exhibit in thick gray regions the Z' mass 
region that accommodates the muon magnetic moment excess. 
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Bounds from Muon Magnetic Moment 


Model 

Overall Sign 

Bound 

U (1) universal 

A a M > 0 

A a M > 0 

gz' = 1, Mz> > 323 GeV 
gz> = 0.5, M z > > 162.5 GeV 

U( 1)b—L 

A a M > 0 

A a M > 0 

gz' = 1, Mz 1 > 323 GeV 
g z > = 0.5, M z > > 162.5 GeV 

^( 1 ) 10+5 

A a M < 0 

A < 0 

gz' = 1, Mz> > 240 GeV 
g z < = 0.5, M z < > 120 GeV 

U( l)d-u 

A a M < 0 

A a M < 0 

gz' = 1, Mz' > 107 GeV 
gz' = 0.5, Mz' > 53.5 GeV 


Table 2. Bounds on the Z' mass rising from the muon magnetic moment for each of the models. 
Those limits are also indicated in Figs. 4-27 as horizontal gray lines. 

5 Collider Bounds: search for new resonances in dilepton events 

The new Z' boson associated with the spontaneously broken U (I)a' gauge group generically 
decays to charged leptons, light jets, top quarks, and, invisibly, to dark matter and/or 
neutrinos. Those decays enable Z' searches with the Large Hadron Collider (LHC) in a 
variety of channels, including dileptons, dijets and rnono-X signatures. As pointed out 
in Ref. [25], the constraints from dijet and monojet searches are competitive only in the 
regime of very weak couplings to charged leptons - e.g. in the “leptophobic” scenario 
studied in [24]. Incidentally, we find that the features of dilepton and dijet events are very 
similar to those presented in [25]. 

Fig. 2 illustrates the decay branching ratios for the Z' into jets (dotted lines), leptons 
(solid lines) and invisible channels (dashed lines). The left-most column assumes M x = 15 
GeV, the middle 500 GeV, and the right-most 1 TeV, while the four rows correspond to 
the four charge structures under consideration: L/(1)universal (first row), U(1)b-l (second 
row), U( l)io+x5 (third row) and U(l)d-xu (fourth row). 

Depending on the U(l)x charge assignment structure, stronger couplings between 
charged leptons and the Z' produce larger branching ratios to charged leptons. This is the 
case e.g. for B — L models, as we can see in the second row of Fig. 2. While a ~ 40% 
branching ratio into charged leptons (including tau leptons decays) is typical for these 
types of models, the other three models considered in this work present smaller branching 
ratios, of order ~ 10%. On the other hand, stronger couplings to quarks lead to larger 
Z' production rates as in the H(l) un i versa i models. The trade-off between production cross 
section and branching fraction to charged leptons determines the relative number of events 
for each model given a specific mass spectrum; the cut efficiencies for dilepton searches are 
expected to be the similar for these models once the Z' mass is fixed. 

Fig. 2 also illustrates that a dilepton search for a Z' is critically sensitive to the model 
details. First of all, for g x = 1, the branching ratio to jets is considerably smaller than 
those found in the leptophobic scenario analyzed in Ref. [24] over the entire Z' mass range 
for M x < 500 GeV, but increases in the Z l mass range where the decay to a DM pair 
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Figure 2. Branching ratios of Z' into jets (dotted lines), leptons (solid lines) and an invisible 
(neutrinos plus DM) mode (dashed lines) as a function of the new gauge boson mass, for DM 
masses M x = 15 GeV (left column), 500 GeV (middle column) and 1 TeV (right column). The first 
row shows the branching ratios of the U(l) un i versa i model, the second row the U(1)b-l, the third 
row the C/(1) 10 +g, and the fourth row the U(l)d- U model. We have fixed gz' = g x = 1 in all plots. 
In all cases there is also a small branching fraction to top quarks, which is not shown. 


is kinetically forbidden. However, constraints from dileptons are much tighter than those 
from dijets even in those cases. This can be understood in terms of the backgrounds 
associated with dileptons and dijets searches for new resonances. While the signal cross 
section for dileptons and dijets from Z' production have similar rates, the backgrounds for 
dijet resonances involve QCD production of jets, which demands much tighter cuts to clean 
up the SM events. On the other hand, the main background for dileptons searches is the 
Drell-Yan process, which is much easier to suppress. 

As in [25], the branching ratio to an invisible final state is larger compared to a lep- 
tophobic scenario due to potential decays to neutrinos. In all models but H(l) un i versa i, the 
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invisible branching ratio can reach the 60% level. This is larger than the typical branching 
ratios of the leptophobic scenario [24]. However, this somewhat larger branching ratio to 
invisible does not make the monojet channel competitive in the case of these models with 
sizable lepton couplings. 

In Ref. [24], Z' masses up to 2.1 TeV for light DM could be ruled out, based on dijet 
and monojet searches at the Tevatron and LHC for g x < 1. Dilepton searches are more 
efficient to exclude regions of parameter space for dark Z' models which present sizable 
lepton couplings. This is the case for the models presented in this work and the model 
independent framework studied in Ref. [25]. 

Dileptons at the 8 TeV LHC - ATLAS search 


To evaluate the constraints from collider searches on the models under consideration, 
we used results from the 8 TeV LHC search for dilepton resonances after 20.3 fb -1 of 
integrated luminosity for the electron sample and 20.5 fb -1 for the muon sample [98]. The 
following Drell-Yan type process was simulated 

pp^Z' ^£~£ + , £ = e,p (5.1) 

plus up to two extra jets using Madgraph5 [99]-FeynRules [100], clustering and hadronizing 
jets with Pythia [101], and simulating detector effects with Delphes3 [102]. Soft and 
collinear jets from QCD radiation generated by Pythia were consistently merged with 
the hard radiation calculated from matrix elements in MLM scheme [103] at appropriate 
matching scales. We adopted the CTEQ 6 L parton distribution functions [104], computed 
at the factorization/renormalization scales pf = Pr = Mz'- 

Signal events were selected according with the same criteria adopted in Ref. [98]: 

pr{e i) > 40 GeV , pri^) > 30 GeV , \rj e \ < 2.47 (5.2) 

Pr( Pi) > 25 GeV , pt{P 2 ) > 25 GeV , \rj,j\ < 2.47 (5.3) 

128 < M u < 4000 GeV (5.4) 

Here £ 1 ( 62 ) is the hardest (second hardest) lepton in the event, and Mu the invariant mass 
of the lepton pair. The signal acceptance times efficiency found in our simulations are 
similar to those presented in [98] and in our previous study [25]. 

All background simulations to pp —>• Z l —>• t~ l + were taken from Ref. [98]. To constrain 
a Z' model we calculated a x 2 statistic (at 95% confidence level) based on Mu measured in 
Ref. [98] in 6 invariant mass bins: 110 — 200 GeV, 200 — 400 GeV, 400 — 800 GeV, 800 — 1200 

GeV, 1200 — 3000 GeV, and 3000 — 4500 GeV. We show in Fig. 3 the number of events 

for the assumed luminosity in each £~£ + invariant mass bin for the total background and 
signal in the f7(l) un i ver sai model (with gz' = g x = 1 and M x = 100 GeV) for various Mz' 
values. The limits we obtain in this work do not take systematic uncertainties into account 
in the fitting procedure, thus our results might be somewhat overestimated. 

The exclusion regions from collider bounds do not vary drastically from one model 
to another in the weak g x coupling regime, as can be seen in Figs. 4-27. This feature 
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Figure 3. The clilepton invariant mass distribution for the total background (dashed line) in six 
bins taken from Ref. [98] and the signal for four different Z' masses in the U( 1) universal couplings 
model: 500 GeV, 1, 2 and 3 TeV. We fixed gz 1 = g x = 1 and M x = 100 GeV. The integrated 
luminosity assumed in this plot is the same as in the ATLAS study [98]. 

is the result of the relative contribution between the Z' production cross section and the 
branching ratio to charged leptons. The branching ratio to light jets, as shown as the 
dotted lines in Fig. 2, illustrates the relative size of the production cross section in each 
model. The CT(1)universal model is expected to produce the largest number of new gauge 
bosons, but precisely because of their stronger couplings to quarks, the branching ratio to 
£+£~ is suppressed, compensating for the larger production cross section compared to the 
other models. For the d — u and 10 + 5 models, the branching ratios to charged leptons 
are larger compared to the universal model, but the couplings to quarks are smaller, again 
compensating each other and rendering the bounds very similar. The B — L model presents 
a much larger branching fraction to leptons which translates into resonance searches probing 
Z' masses roughly ~ 300 GeV larger than the other models under consideration. 

The other prominent feature we observe in Figs. 4-27 is how the bounds become in¬ 
creasingly weaker as g x increases. As seen in Table 1, in d — u models neutrinos do not 
couple to the new force. As a result, the branching ratio to invisible states arises solely 
from decays to dark matter, thus enhancing the sensitivity of this model to the precise 
value of the coupling g x . The 10 + 5 model also introduces rather weak couplings between 
the Z l and neutrinos compared to the universal and B — L models. This explains why the 
bounds on Mz> are nearly constant as a function of g x in the universal model, but weaken 
in a much more pronounced manner for the d — u model. However, in the small g x regime, 
all these models are excluded for Z' masses below ~ 2.6 TeV and dark matter masses from 
8 GeV to 5 TeV, except for the universal model where masses below ~ 3 TeV are excluded 
for the same dark matter masses. Our results for g x ~ 0 compare favorably to the ATLAS 
bounds on other Z' models [98]. For larger couplings, the bounds are softened, as discussed 
above, reaching ~ 1.8 TeV in the d — u models with g x = 1 and M x = 8 GeV and ~ 2.4 
TeV for the same couplings with M x = 1 TeV. We point out that the limits from LEPII 
derived in [ 68 ] are not applicable here due to the existence of invisible state. 
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Collider Limits 


Model 

Mass 

Coupling 

Bound 

U ( 1 ) universal 

8 GeV- 1 TeV 

8 GeV- 1 TeV 

9z' = 1 
9z' = 0.5 

g x < 0.1, M z > > 2570 GeV; g x ~ 1, M z > > 2400 GeV 
g x < 0.1, M z < > 2196 GeV; g x ~ 1, M z > > 2100 GeV 

^( 1 ) 10+5 

8 GeV - 1 TeV 

8GeV - 1 TeV 

9z> = 1 
9Z' = 0.5 

g x < 0.1, M z > > 2620 G eV;g x ~ 1, M z > > 2400 GeV 
g x < 0.1, M z > > 2480 G eV;g x ~ 1, M z > > 2040 GeV 

U( 1)b-l 

8 GeV - 1 TeV 

8 GeV - 1 TeV 

9z> = 1 
9Z' = 0.5 

g x < 0.1, M z > > 3000 GeV ;g x ~ 1, M z > > 2930 GeV 
g x < 0.1, M z > > 2570 GeV;g x ~ 1, M z > > 2280 GeV 

U{ l) d - u 

8 GeV-1 TeV 

8 GeV-1 TeV 

9z' = 1 
9Z' = 0.5 

g x < 0.1, M z , > 2640 G eV-g x ~ 1, M z , > 2300 GeV 
g x < 0.1, M Z i > 2430 GeV;g x ~ 1, M z > > 1800 GeV 


Table 3. Collider limits on the Z' masses for the models in study. These limits are also reproduced 
in the plots of Figs. 4-27. We point out that the bounds for g x ~ 1 are actually sensitive to the 
DM mass, we thus quote the most conservative ones. 


We also should comment that we are ignoring interference effects with the SM Z boson 
and photon. For the most part of the parameters spaces of the models under study, the 
zp width is sufficiently narrow to justify that. Yet, even for larger widths, these effects 
are expected to be small for Z' masses close to the bounds that we found. Moreover, 
ignoring the interferences gives us more conservative bounds, once the interference among 
the neutral vector states is constructive. 

6 Results 

We have outlined a comprehensive study of four different U(l)x realizations, in the context 
of new physics contributions to the muon magnetic moment, collider searches, and direct 
and indirect dark matter searches. In addition to considering different, theoretically well 
motivated fermion charge assignments, we focused on two values for the gz' coupling, 1 
and 0.5. The muon magnetic moment and collider limits are summarized, respectively, in 
Tables 2 and 3. The collider bounds were derived using dilepton search results from ATLAS, 
Ref. [98]. We emphasize that the limits from LEPII derived in [ 68 ] are not applicable here 
due to the existence of dark matter. 

As for the muon magnetic moment bounds, we utilized the public code described 
in [59]. For the direct and indirect detection limits we compared our analytical results 
with micromegas [105] concerning the thermal relic abundance calculation and utilized the 
PPPC4DMID code [83] to calculate the gamma-ray spectra. 

In this section we present a summary of our results, combining all constraints outlined 
above, for a variety of dark matter particle masses, and for all the charge assignments 
described, on the plane defined by the Z' mass and the DM-Z' coupling g x . Each plot 
focuses on a particular value for gz 1 = 0.5, 1 and M x = 8 , 15, 50, 500, 1000 and 5000 
GeV, and a selected charge assignment structure. The broad range of DM masses allows 
us to highlight the importance and complementarity across different search strategies. 

General comment # 1: We employ the same color scheme for all figures: The red 
(purple) regions indicate the parameter space ruled out by LUX SI (XENON SD) direct 
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Figure 4. Results for the 17(1) universal model. Left: m x = 8 GeV, g Z ' = 1; Right: rn x = 8 GeV, 
gz' = 0.5. The blue horizontal line indicates the ATLAS LHC bound: everything below the curve 
is ruled out. The gray horizontal line is the lcr bound from the muon magnetic moment by forcing 
the correction to lie within the error bar. The thick gray band delimits the Z' that accommodates 
the muon anomalous magnetic moment excess. The red (green) regions are ruled out by the LUX 
spin-independent direct detection results, while the green region is ruled out by Fermi-LAT dwarfs 
PASS8 limits. The black curve indicates the region of parameter space that features a DM relic 
density matching the universal DM abundance. The cyan shaded region corresponds to the violation 
of the perturbative limit, T z > > M Z '/2 . In the left graph (g Z ' = 1), the whole parameter space 
violates this pertubative limit, therefore the dilepton limit does not apply since it is based on the 
narrow width approach. 




Figure 5. Same as in Fig. 6, for the 17(1)universal model. Left: m x = 15 GeV, g Z ' = 1; Right: 
m x = 15 GeV, g Z ' = 0.5. 


dark matter searches (Sec. 3.1). The black curve indicates the region of the parameter 
space where the thermal relic density of % matches the observed dark matter abundance, 
while the green shaded region is ruled out by gamma-ray observations of Dwarf satellites 
(Sec. 3.2). The horizontal gray line represents the muon magnetic moment limit (Sec. 4), 
whereas the gray band delimits the region of parameter space that accommodates the 
g-2 excess. In all cases, the region compatible with a new physics interpretation of the 
anomalous magnetic dipole moment of the muon is ruled out by other other searches. 

General comment # 2: The blue curves indicate the minimal M z > compatible with 
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Figure 6. Same as in Fig. 6, for the 17(1)universal model. Left: rn x = 50 GeV, gz' = 1; Right: 
m x = 50 GeV, gz' = 0.5. 




Figure 7. Same as in Fig. 6, for the 17(l) un iversai model. Left: m x = 500 GeV, gz 1 = 1; Right: 
m x = 500 GeV, gz 1 = 0.5. 




Figure 8 . Same as in Fig. 6, for the 17(1)universal model. Left: rn x = 1 TeV, gz> = 1; Right: 
m x = 1 TeV, gz 1 = 0.5. 
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Figure 9. Same as in Fig. 6, for the U(l) universal model. Left: 5 TeV WIMP, gz' = 1; Right: 5 
TeV WIMP, gz> = 0.5. Notice that a DM mass above a few TeV is required to be consistent with 
current limits and reproduce flh 2 =0.1 

null ATLAS searches in the dilepton channel (Sec. 5). The cyan region represent the 
parameter space where the total width of the Z' exceeds the perturbative limit, Tz> > 
Mz< / 2 , invaliding the dilepton bounds which are based on the narrow width approximation. 
However, one should notice that the direct detection limits are most relevant in this regime 
leaving our general conclusions intact. In particular, this perturbative bound on the total 
width of the Z' is exceeded for all cases where we consider the C/ (1) un i versa i model with 
gz 1 = 1. However, for completeness, we choose to still display these results. 

General comment # 3: Concerning the collider bounds, one can see from Table 1, 
that some gauge structures induce stronger Z'-leptons couplings, and consequently larger 
branching ratios to charged leptons, such as the B — L model. Moreover, gauge structures 
such as t/ (1)universal lead to sizable Z'-quarks couplings and thus large production rates. 
The balance between production rate and branching ratio to charged leptons sets the 
number of dilepton events at each invariant mass bin, and therefore our limits. 

General comment # f: We notice that in the absence of vector couplings to quarks, 
RG running from the Z' mass scale down to nucleon energy scales effectively induces spin- 
independent interactions, as pointed out and investigated in Ref. [106]. The resulting 
limits, which we do not compute here, might be comparable or even stronger than those 
from spin-dependent searches. 

We now summarize the results presented in Figs.4-27: 

• H ( 1 )universal • I n Figs. 4-9 we show the exclusion region for several dark matter masses 
and conclude that only for DM masses larger than ~ 1 TeV can one obtain models 
with the right relic abundance which evade both direct detection and collider con¬ 
straints. For gz' = 1(0.5), collider bounds roughly require Mz> > 2570(2200) GeV. 
This excludes this setup from providing an explanation to the ATLAS diboson ex¬ 
cess [60], since a large gz 1 is needed to produce a large enough number of events 
[62] . Note that these limits are marginally sensitive to the magnitude of the Z'-dark 
matter coupling (g x )- When g x becomes sufficiently large the limits weaken because 
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the DM branching fraction increases and the charged lepton one decreases. Keeping 
gz' = 1, in particular we find Mz> > 2400 GeV, g x ~ 1. A summary of our limits is 
provided in Table 3. 

In this model the correction to g-2 has the right sign, however the region of parameter 
space that accommodated the g-2 excess is ruled out as one can see. For g x ~ 10 -2 — 
10 -1 the direct detection limits the most constraining ones, whereas the indirect 
detection bounds are always weaker. In summary, DM masses above few TeV are 
required to match the universal DM abundance and obey the existing bounds. 

• C/(1)b-l : The combination of g-2, collider, direct a indirect dark matter detection 
limits are depicted in Figs. 10-15. The collider limits result in Myj > 3 TeV for 
gz> = 1 and Mz' > 2570 GeV for gz> = 0.5. Again, this is incompatible with the 
ATLAS diboson excess. See Table 3 for different coupling choices. Similarly to the 
1 /(1) universal model, in principle it could accommodate the g-2 excess, but the region 
of parameter space is excluded by collider data, and DM masses larger than few TeV 
are needed to reproduce Qh 2 = 0.1 while obeying the current limits. 

• t/(l) 10+ g: We present our findings for this gauge model in Figs. 16-21. Keeping 
g x ~ 0.1, the collider bounds are in the ballpark of Mz> <1 2.6 TeV for gz> = 1 and 
Mz> > 2480 GeV for gz> = 0.5. For g x ~ 1 those limits drop to Mz> > 2400 GeV 
and Mz< > 2040 GeV respectively. In this model the usual leading contribution 
to the SI cross section vanishes at tree level, since the vector couplings to valence 
quarks vanishes. However, spin-dependent limits from XENON100 are still relevant, 
and are complementary to the Fermi Dwarf limits for certain DM masses, although 
generically weaker. As mentioned in comment # 4 above, loop effects might produce 
sizable spin-independent scattering as well. It is clear from Fig. 19 that one can have 
a Dirac dark matter WIMP with mass of 500 GeV or larger in this model. 

• U(l)d- U - Our results are summarized in Figs. 22-27 and Table 3. The importance 
of taking a dark matter complementarity approach is clear in this model, since for 
g x < 1, the dilepton limits are the strongest, whereas for g x > 1, the direct and 
indirect detection bounds are the leading ones. The dilepton limits are again very 
stringent, excluding Z' masses below 2640(2430) GeV for gz’ = 1(0.5). For gz' 
of order of unit those limits drop to 2430(1800) GeV respectively. In this model, 
similarly to the previous case, for masses of 500 GeV one can accommodate a Dirac 
fermion as DM 1 . 

7 Diboson Excess 

We show in this section that it is possible to accommodate the recent ATLAS diboson 
excesses [60] within the U(l)d- U model discussed above. The diboson excesses in the 

2 We point out our results only apply for Dirac fermions which were thermally produced in the early 
universe. For recent discussion concerning non-thermal DM production see Ref. [1, 107-114]. 
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WZ,WW and ZZ channels reported by the ATLAS Collaboration are well fitted 
by resonances whose peaks are in the ballpark of 2 TeV, potentially implying the 
detection of a new particle. Since the tagging selections for each mode used in the 
analysis are relatively poor at this stage (~ 20%) it is hard to conclude that one 
resonance is responsible for the excesses in all channels. Indeed, it is possible that a 
single 2-TeV particle such as a Z l might account for only one part of the channels, and 
that the peaks in the other channels are contaminations due to incomplete tagging 
selections, see e.g. [62], 

In the context of U(l)d- U models, we find that we can accommodate the excess while 
reproducing the right thermal relic abundance for the Dirac fermion dark matter 
candidate we discuss here, avoiding at the same time dilepton, direct and indirect 
detection bounds. We incorporated a Z'WW interaction in the U{l)d~ u setup, pro¬ 
portional to gz'8, where 8 is the Z — Z' mixing angle, bounded to be smaller than 
10 -3 [68] due to constraints from SM Z properties. We then computed the relic abun¬ 
dance, dilepton bounds and the pp —> Z' —> W W production cross section, finding 
that for Mz' ~ 2 TeV, gz r = 1; g x ~ 1, m x = 500 GeV; 0 = 10 -3 , we can reproduce 
the right abundance ( VLh 2 ~ 0.1) and obtain a(pp —> Z' — > WW) ~ 30 fb, which 
accommodates the excess. 

We note that we employed m x = 500 GeV because for lighter masses we cannot 
simultaneously get the right abundance and avoid dileptons exclusion limits. Fur¬ 
thermore, for m x > Mz* the branching ratio into WW is sizable, inducing a WW 
production cross section too large to reproduce the excess. However, we find that 
it is possible to accommodate heavier dark matter, with masses up to ~ 1 TeV by 
adjusting the production cross section of a 2 TeV Z’ in order to explain the diboson 
excess, and at the same time we respect all the other bounds considered in this work, 
since the inclusion of a WW decay channel weakens the dilepton collider bounds, 
thus opening a broader region of the allowed parameters space which could explain 
the ATLAS excess. 

8 Conclusions 

In this work we studied several anomaly-free 1/(1) gauge structures that arise in several, 
well-motivated extensions of the standard model, namely JJ{1)b-l, U(l)d- U , U(l) urL i versa i, 
and 1/(1) 10+ g. By postulating a Dirac fermion x coupled to the Z' as the dark matter 
particle candidate, we computed the thermal relic abundance, annihilation cross section in 
the low velocity limit, and spin-dependent and -independent scattering cross sections off of 
nuclei. We then used the current bounds from dwarfs observations using Fermi-LAT data 
and from null results reported by LUX and XENON to set constraints on the Mz r versus 
g x parameter space, where g x is the Z'-DM coupling. 

We derived dilepton limits from null searches with ATLAS, which provide the strongest 
constraints for all those models for g x < 10” 2 , with direct detection limits leading for larger 
couplings. Dilepton constraints rule out Z l masses up to ~ 3 TeV. In general, the collider 
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limits depend quite sensitively on the assumed gauge charge structure, but not as much 
on the g x coupling. Albeit, the collider limits for the U(l)d- U gauge group do weaken 
significantly for g x > 1. 

In no case we find large enough contributions to the anomalous magnetic moment of the 
muon large enough to explain the observed anomaly compatible with collider constraints, 
typically the combination of collider and direct detection limits forces the mass of the dark 
matter particle to values larger than a TeV, assuming thermal production in the early 
universe, unless very suppressed couplings to leptons are in place (as is the case for the 
C/(l)io +5 model, for which the limits is around 0.5 TeV). Relaxing the thermal production 
requirement, but enforcing a thermal abundance at most as large as the 2 a upper limit on 
the inferred universal DM density does not change this conclusion. Underabundant models 
generically prefer small values for Mz 1 and large values for g x , but are severely constrained 
by indirect detection constraints. 

Finally, we addressed the question of whether the models under investigation could 
accommodate the recently reported ATLAS diboson excess; we found that for large enough 
Z' couplings a Z' mass on the order of 2 TeV is generically excluded for all but the U(l)d- U 
model by a combination of LHC and direct detection constraints. For the U{l)d-u nrodel 
a 2 TeV Z’ and a thermal relic Dirac fermion dark matter particle with a mass of 500 GeV 
appears to be a possibility, although further investigation is needed to probe whether the 
required number of events could be produced in this scenario. 
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